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Abstract 

We develop a theory to represent dislocated single crystals at the meso- 
scopic scale by considering concentrated effects, governed by the distribu- 
tion theory combined with multiple- valued kinematic fields. Our approach 
gives a new understanding of the continuum theory of defects as devel- 
oped by Kroner (1980) and other authors. Fundamental 2D identities 
relating the incompatibility tensor to the Frank and Burgers vectors are 
proved under global strain assumptions relying on the geometric measure 
theory, thereby giving rise to rigorous homogenisation from mesoscopic to 
macroscopic scale. 

1 Introduction 

Dislocations can be considered as the most undesirable and resistant class of defects 
for several kinds of single crystals (Maroudas and Brown, 1991; Jordan et al., 2000). 
Therefore, in order to improve crystal quality, the development of a relevant and ac- 
curate physical model represents a key issue with a view to reducing the dislocation 
density in the crystal by acting in an appropriate way on the temperature field and the 
solid-liquid interface shape during the growth process (Van den Bogaert and Dupret, 
1997). 

However the dislocation models available in the literature, such as the model of Alexan- 
der and Haasen (1968, 1986), are often based on a rather crude extension of models 
previously developed for polycrystals (such as usual metals and ceramics). In this 
case, some particular features of single crystals, such as material anisotropy or the 
existence of preferential glide planes, can be taken into account up to some extent, but 
the fundamental physics of dislocations in single crystals cannot be captured. In fact, 
dislocations are lines that either form loops, or end at the single crystal boundary, or 
join together at some locations, while each dislocation segment has a constant Burgers 
vector which exhibits additive properties at dislocation junctions. These properties 
play a fundamental role in the modelling of line defects in single crystals and induce 
key conservation laws at the macro-scale (typically defined by the crystal diameter). 
On the contrary, no dislocation conservation law exists at the macro-scale for poly- 
crystals since dislocations can abruptly end at grain boundaries inside the medium 
without any conservation law holding across these interfaces. 
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Aware of these principles and of the pioneer works of Voherra (1907) and Cosserat 
(1909), Burgers (1939), Eshelby, Frank and Nabarro (1951, 1956), Kondo (1952), Nye 
(1953), and Kroner (1980) among other authors (Bilby, 1960; Mura, 1987) consider a 
tensorial density to model dislocations in single crystals at the macro-scale, in order to 
take into account both the dislocation orientation and the associated Burgers vector 
(cf the survey contributions of Kroner (1980, 1990) and Kleinert (1989)). However, 
in these works, the relationship between macro-scale crystal properties and the basic 
physics governing the nano-scale (defined by the inter-atomic distance) is not com- 
pletely justified from a mathematical viewpoint. Therefore, to well define the concept 
of tensorial dislocation density, we here introduce an additional scale to the macro-and 
nano-scales, viz. the meso-scale as defined by the average distance between disloca- 
tions. The laws governing dislocation behaviour are modelled at the nano-scale, while 
the meso-scale (defined from the nano-scale by ensemble averaging or by averaging over 
a representative volume (Kroner, 2001)) defines the "dislocated continuous medium", 
where each dislocation is viewed as a line and the interactions between dislocations 
can be modelled while the laws of linear elasticity govern the adjacent medium. 
The present paper focuses on meso-scale modelling with a view to clarifying the ho- 
mogenisation process from meso- to macro-scale. Since dislocations are lines at the 
meso-scale, concentrated effects, as governed by the distribution theory (a key ref- 
erence is here Schwartz (1957)), must be introduced in the mesoscopic model. In 
addition, since integration around the dislocations generates a multiple-valued dis- 
placement field with the dislocations as branching lines, multivalued functions must 
be considered (Almgren, 1986; Knopp, 1996; Remmert, 1996). This combination of 
distributional effects and multivaluedness is a key feature of the dislocation theory 
at the meso-scale but, unfortunately, the difficulties resulting from this mathematical 
association have not well been addressed so far in the literature (Thom, 1980). As an 
example, non-commuting differentiation operators are freely introduced without any 
justification by Kleinert (1989). Therefore, the principal objective of this paper is to 
provide a strong mathematical foundation to the meso-scale theory of dislocations, 
showing how the distribution and geometric measure theories can be correctly used 
with multiple-valued fields. In particular, the application limits of Stokes' theorem 
are clarified. For the sake of generality, disclinations, which represent a second but 
rarer kind of line defect, with in addition a multiple-valued rotation field, are here 
considered together with dislocations. 

After homogenisation from meso- to macro-scale, no concentrated effects remain any- 
more present in the macroscopic model, which consists of a set of evolution PDE's 
governing scalar or tensorial defect density fields in the framework of elasto- or visco- 
plasticity (Kratochvil and Dillon, 1969). However, it should be pointed out that 
homogenisation from meso- to macro-scale has no meaning for multiple-valued fields 
such as displacement and rotation, since this operation is exclusively allowed for addi- 
tive (or extensive) fields such as stress, energy density or heat flux. This consideration 
becomes obvious when homogenisation is defined by an ensemble averaging procedure, 
since multiple-valued fields are mathematically defined as extended functions which 
cannot be added since their "domains" depend on the defect line locations. This issue 
is clarified in the present paper. Moreover, since the macroscopic displacement and 
rotation fields are not defined as ensemble averages of their mesoscopic counterparts, 
no unique privileged reference configuration can be defined at the macro-scale for sin- 
gle crystals with dislocations. Having in mind that displacement and rotation fields 
are defined with respect to the selected reference configuration (which can be, or not, 
defect free), the invariance laws governing the behaviour of single crystals with line 
defects at the macro-scale are constructed in accordance with this observation. 
In the literature the macroscopic dislocation density is classically defined as the curl 
of the plastic distortion (Head et al., 1993; Cermelli and Gurtin, 2001; Gurtin, 2002; 
Koslowski et al., 2002; Ariza and Ortiz, 2005), following a postulated distortion decom- 
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position into elastic and plastic parts. However, this approach cannot be rigorously 
justified (contrarily to the strain decomposition) since elastic and plastic rotations 
cannot be set apart. In contrast, the present paper introduces the macroscopic dis- 
location density by homogenisation of well-defined mesoscopic fields, under precise 
geometric-measure model assumptions, from which the distortion decomposition is 
obtained together with its relationship with the dislocation density. Since dislocations 
and disclinations represent body torsion and curvature, respectively, these concepts 
also appear as macroscopic counterparts of well-defined mesoscopic defect measures. 
The present paper is restricted to the 2D theory. Extension to the 3D case is under 
investigation. A complete link between the mesoscopic and macroscopic behaviours 
of single crystals with line defects should be derived from these developments. In 
Section 2, the scaling analysis summarised in this introduction is detailed and the 
basic concepts used to represent the dislocated continuous medium are introduced. 
Classical invariance theorems are recalled in Section 3. In Section 4, the 2D distribu- 
tional theory of the dislocated continuous medium is established in the case of a single 
dislocation, while Section 5 treats the more general case of an ensemble of isolated 
dislocations. Finally, Section 6 introduces the non-Riemannian macroscopic body by 
homogenisation from the mesoscale, and conclusions are drawn in Section 7. 

2 Multiscale analysis of dislocations 

To address the modelling of elastic single crystals with line defects, the various scales 
relevant for matter description and their interrelations are briefiy reviewed. 

• At the nano-scale the characteristic length is the interatomic distance. At time 
t, the body is referred to as TZ**{t) and the reference body is a perfect lattice 
TZ*o*. 

• At the meso-scale the characteristic length is the average distance between two 
neighbour dislocation lines. At time t, the body is referred to as TZ*{t), to be 
interpreted as a random sample corresponding to a given growth experiment. 
The reference body TZq is a perfect crystal, i.e. a body without dislocations or 
disclinations. 

• At the macro-scale the characteristic length is the diameter of the crystal and 
the body TZ{t) has a physical meaning related to TV it) and 7l**{t) in terms of 
ensemble average; the reference body TZo can be, or not, a perfect crystal. 

2.1 Nano-scale analysis: crystalline lattice 

Given a dislocation in the general sense, the atomic arrangement TZ**{t) differs from 
the reference lattice TZq* , but however the atom displacements are not uniquely defined 
(Kleinert, 1989). More exactly, a discrete multivalued mapping Xi := Xi*{^i) where 
i = 1,2 or 3, is defined with Xi £ TZq* and Xi £ TZ**{t). In general, the dislocation 
position cannot be determined precisely at the atomic level since several dislocation 
locations in the actual crystal can be associated with the same picture of the atom 
positions. In fact the defect should be understood as located inside a nanoscopic 
lattice region. Let us insist on the fact that there is no way to uniquely define the 
displacement field. Indeed any atom of TZq* can in principle be selected to define 
the displacement of a given atom of TZ**{t) which is therefore a multivalued discrete 
mapping. This remark also makes sense at higher scales. 

2.2 Meso-scale analysis: dislocated continuous medium 

This scale is the one on which this paper focuses, in the framework of 2D linear 
elasticity. Let us here describe some general and basic field properties at the meso- 
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scale level. Details are given in the forthcoming sections. 

• The displacement field is a multivalued function such that for any point Xi G TZq 
one has 

u*[Xi) = Xi-Xi, with Xi--x*i{Xi)^ 

and where Xa^Xi) is a multivalued mapping from TZ*-, to 71* (t). Displacement 
multivaluedness represents an important difficulty to address at the meso-scale 
in dislocation modelling. As opposed to multiple-valued fields, single-valued 
fields will also be called uniform. 

• The strain will be denoted by £*j. In general, the Lagrange deformation tensor 
is given by 



i 



^iJ •= 2^9ju* + diu* + djul^diu^), with djU* := 



with the classical indicial notations and summation convention used. In the 
sequel, linear elasticity will be assumed and hence the nonlinear terms are not 
taken into account at the meso-scale. This fundamental hypothesis relies on the 
assumption^ that all nonlinear deformation effects take place around the disloca- 
tion in a nano-scale region whose diameter is small compared to the meso-scale 
characteristic distance. Therefore, using a singular perturbation asymptotic 
treatment, the nonlinear effects become concentrated inside the defect line at 
the meso-scale and hence the strain can be assumed to be the single-valued 
linear symmetric tensor given by 



^ij '■= T^idjUi + diu*) 



outside the defect line and arbitrarily set to on the defect line (noting that 
concentrated deformation effects inside this line do not play any role in displace- 
ment integration at the meso-scale). 

• The infinitesimal rotation tensor is a possibly multiple-valued field given by 
LOij :— ^{djU* — diU*) with the associated rotation vector given by 

1 . 1 „ * 

and the identity uj*j = —eij^u^. The Frank and Burgers vectors Q,*^. and B* 
associated with a defect line are commonly defined as functions of the jumps 
of ujI. and u* around this line. From Weingarten's theorems, these vectors are 
shown as invariants of the defect line. 

The following geometric tensors are also introduced: 
Definition 2.1 

DISCLINATION DENSITY: 9 



DISLOCATION DENSITY: A. 



DISPLACEMENT JUMP DENSITY: 



CONTORTION: K*^ — a*^ - -a^^^Sij , (2.4) 



^n*s,L, (2.1) 

= B*S,L, (2.2) 

= A*j -I- eji^Q*i{xm — xom); (2.3) 
1 



where xom ts a reference point for rotation and displacement integration. 



^In practise this assumption is certainly valid in single crystal growth. 



4 



Here, the symbol SiL is used to represent the concentrated vectorial measure density 
on the defect region L. In particular, when L is a rectifiable curve, SiL is equal to TiSL 
with Ti the unit tangent vector to L and Sl the one-dimensional Hausdorff measure 
density concentrated on L. 

The disclination and dislocation density tensors Q*j and A*j are measure densities (cf 
Evans and Gariepy, 1992; Mattila, 1995) related to the so-called strain incompatibility 
riij to be defined later. In general at the meso-scale a dislocation or a disclination 
is a defect line (i.e. a singular line for the strain) to which non-vanishing Burgers 
and/or Frank vectors are attached. The tensors A*^ and Q*j are basic physical tools 
to model defect density at the meso-scale while 77* plays a key role to understand 
their behaviour. The displacement jump density and mesoscopic contortion (or lattice 
curvature) tensors a*j and combinations of these basic density tensors, with 

a*j = A*j when the disclination density tensor vanishes. 

2.3 Macro-scale analysis: continuous medium 

At this level, a point Xi of the actual body TZ{t), Xi = Xi{-^i) where X £ TZo will be 
called a material point to be understood as a certain volume of matter of mesoscopic 
size located around the point Xi. In order to define macroscopic concepts such as 
temperature or stress, one needs to give a meaning to the temperature and stress at 
any point. The rigorous definition is obtained from an ergodicity argument and hence, 
at the macroscopic level the fields on TZ{t) are defined as ensemble averages of the 
fields defined on TZ*{t). By this operation these fields are smoothed, which means that 
concentration effects at the meso-scale level along the defect lines are erased. To this 
end, a weak limiting procedure (or homogenisation) is needed in order to define the 
dislocation and disclination densities Aij and Oij at the macro-scale level from the 
knowledge of the meso-scale fields A*j and Q*j. 

Remark 2.1 In this context, the reference body TZo is basically arbitrary and can, 
or not, be a perfect crystal. Indeed, at the macro-scale, the displacement Ui must be 
a single-valued function, while the displacement field u* is multivalued at the meso- 
scale. Therefore, u* cannot be considered as belonging to a (linear) Banach space 
(single-valued functions can be added since they share the same domain, whereas a 
multivalued function is defined on its specific Riemann foliation and cannot be added 
to a multivalued function defined on another Riemann foliation). Consequently the 
ensemble averaging procedure is forbidden for multivalued fields such as u* and hence 
Ui is not the ensemble average ofu*. It should also be observed that removing the field 
multivaluedness by performing appropriate cuts is of no use here, since by derivation 
these cuts introduce arbitrary distributional contributions without physical meaning. In 
general, it is important to make it clear that the only fields which can be obtained at 
the macro-scale by ensemble averaging from the meso-scale are the so-called extensive 
fields associated with additive physical properties (such as specific mass, stress, specific 
internal energy... and the dislocation and disclination densities). 

3 Multiple-valued fields and line invariants; dis- 
tributions as a modelling tool at the meso- 
scale 

Notations 3.1 In the following sections, the assumed open domain is denoted by Q 
(in practise but not necessarily ^1 is bounded), the defect line(s) are indicated by L £ Q., 
and Ql is the chosen symbol for ^l\L, which is also assumed to be open without loss of 
generality. Focussing on the meso-scale, symbol x or Xi denotes a generic point of the 
defect line(s), x or Xi is a generic point of and xa or xoi is a given fixed reference 
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point of Q,L- When x and x are used together, x denotes the projection of x onto the 
defect line(s) L in a appropriate sense and Vi := Vi{x,x) is the unit vector joining x 
to X. The symbol 0e is intended for a set of diameter 2e enclosing the region L. More 
precisely, 0e is defined as the intersection with Q of the union of all closed spheres of 
radius e centred on L: 

0, := y B[x,e]. 

In case L is an isolated line, 0^ is a tube of radius e enclosing L. In the sequel, 
considering a surface S of crossed by L at x and bounded by the curve C , symbols 
dC , dL, and dS will denote the ID Hausdorff measures on C and L, and the 2D 
Hausdorff measure on S, respectively, with aj and tj the unit tangent vectors to C 
at X and to L at x (when they exist). In some cases (having fractal curves in mind) 
the symbols dxk and dSi ~ eijkdx''P dx''^'^ will stand for infinitesimal vectors oriented 
along C and normal to S, respectively, with m addition dCiix) := eimndxmTn denoting 
an infinitesimal vector normal to C when Tn — r„(a;) exists. 

Assumption 3.1 (Mesoscopic elastic strain) Henceforth we will assume that the 
linear strain £^„ is a given symmetric C°°(flL, R'^^'^) -tensor prolonged by on the line 
L, -integrable on SI and compatible on SIl. In other words, the equality 

eqlrnekpndldpS^n = (3-1) 

is assumed everywhere on^L. 

3.1 Distributional analysis of the multiple- valued fields 

In general, a multivalued function from Ql to consists of a pair of single- valued 
mappings with appropriate properties: 

F ^Ql and F ^ R'^ , 

where F is the associated Riemann foliation (Almgren, 1986; Knopp, 1996; Remmert, 
1996). In the present case of meso-scale elasticity, we limit ourselves to multivalued 
functions obtained by recursive line integration of single-valued mappings defined on 
flL. Reducing these multiple line integrals to simple line integrals, the Riemann folia- 
tion shows to be the set of equivalence path classes in Ql from a given xo £ with 
the homotopy as equivalence relationship. Accordingly, a multivalued function will be 
called of index n on Ql if its n-th differential is single- valued on SIl. No other kinds 
of multifunctions are considered in this work, whether L is a single line or belongs to 
a more complex set of defect lines (with possible branchings, etc.). 

Notations 3.2 The notation 9^"' is used for partial derivation of a single- or multiple- 
valued function whose domain is restricted to Hl- Locally around x G Ql, for smooth 
functions, the meanings of dj"^ and the classical dj are the same, whereas on the 
entire Q, the partial derivation operator dj only applies to single-valued fields and must 
be understood in the distributive sense. A defect-free subset U ofQis an open set such 
that [/ n 1/ = 0, in such a way that dj"'' and dj coincide on U for every single- or 
multiple-valued function of index 1. 

In the following essential definition the strain is considered as a distribution on Q: 

Definition 3.1 [Frank tensor] The Frank tensor dm^^k is defined on the entire domain 
r2 as the following distribution: 

dmiJ-'k •= ^kpqdpEqm, (3.2) 
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in such a way that 



< dmi^l.,(fi >:= - ekpqSqmdpipdV, (3.3) 

in 



with (fi a smooth test-function with compact support in Q. 

In fact, the tensorial distribution dmcol is the finite part of an integral when acting 
against test-functions. Indeed, since dp£q^ might be non-L^ (r2)-integrable in view 
of its possibly too strong singularity near the defect line, instead of being directly 
calculated as an integral, < €kpqdp£q„^, if > must be calculated on S7 as the limit 

lim / CkpqdpEq^^pdV + / ekpq£gmV>dSp ) , (3.4) 

where the second term inside the parenthesis is precisely added in order to achieve 
convergence. One readily sees after integration by parts that expression (3.4) is equal 
to Eq. (3.3) provided lim f2 \ 0e = Ql (which is a general hypothesis limiting the ac- 

ceptable defect lines together with the assumption that L is of vanishing 2D Hausdorff 
measure). Considering the possibly multivalued (with index 1) rotation vector ujI, it 
should be observed from Definition 3.1 that 9m w J — dm^ujl on Ol. This results from 
the classical relationship provided by elasticity theory between infinitesimal rotation 
and deformation derivatives. However, dmUjl is defined by Eq. (3.2) as a distribution 
on f2 and therefore concentrated effects on L and its infinitesimal vicinity are added 
to dm^Lol, justifying the use of the symbol dmUoX instead of dmi-^k without giving to 
dm the meaning of a derivation operator. In particular, it may be observed that the 
identical vanishing of dm^Lol on Ql does not necessarily imply that the distribution 
dmi^k vanishes as well. In fact from Eq. (3.4), it can be shown in this particular case 
that 

< dml^l,!^ >= Vaa I CkpqSqm'pdSp ^ - ekpqSqmdpipdV, (3.5) 

which is generally non- vanishing. Finally, as soon as the definition of the tensor 
distribution dmiJk is given, so are the distributional derivatives of dmi^l- 

< dldmOJkjf >= - < dmi^t,dllf >= / ekpn£mndpdltpdV. (3.6) 



n 



3.2 Rotation and displacement vectors 

The rotation vector is defined from the knowledge of the linear strain together with 
the rotation at a given reference point xq. From this construction follows an invariance 
property of as a multifunction (recalling that multivaluedness takes its origin from 
the existence of a defect line which renders the strain incompatible on the entire Q,). 
Starting from the distributive Definition 3.1 of dmi^lj the differential form dm^^kd£,m. 
is integrated along a regular parametric curve F C ^II with endpoints xo,x £ For 



r 



where #r is the equivalence class of all regular curves homotopic to F in SIl. Indeed, 
from strain compatibility in Ol, i.e. from relation (3.1), it is clear that is a function 
of #F only. Consider now a regular parametric loop C (in case C is a planar loop, it is 
called Jordan curve) and the equivalence class #C of all regular loops homotopic to C 
in Ol. Here, the extremity points play no role anymore and two loops are equivalent 
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if and only if they can be continuously transformed into each other in JIl. The jump 
of the rotation vector ivl along #C depends on #C only and is defined as^ 

H] = H]i#C) = / dmiuU^n.. (3.7) 
JC 

The following developments address the displacement field multivaluedness as a mere 
consequence of strain incompatibility. The procedure defining the displacement vector 
from the rotation vector by means of line integrals is classical in linear elasticity. The 
following tensor plays in the construction of the displacement field a role analogous to 
dmUil in the construction of the rotation field: 

Definition 3.2 [Burgers tensor] For a selected reference point xo £ Hl, the Burgers 
tensor is defined on the entire domain f2 as the distribution 

dibl{x;xo) ■- £ki{x) + Ckpqixp ~ xop)diuj*{x). (3.8) 

The Burgers tensor can be integrated in the same way as the Frank tensor along any 
parametric curve F, providing for selected xo, ujQk and Mq* the multivalued displace- 
ment vector u*). of index 2: 

ul=ul{ij^T,U)l,ul) = ul^+ eklrn'^*l{x\'^){Xrn - X0„t) + j dlbl{^)d^l, 

which is a function of #F only. It may be observed that dib^. and the vector b^ = 
ul — ekimiOi{xm — xom) are related in the same way as dmUil and wj, including the 
fact that dib^ = ^'^'fej on Ql- The jumps of b^ along ^^^C and of u^. at x along #C 
(which depends on #C only) are defined as 

[bl]{#C;xo) = [ul]{x;#C;xo) - ekim[i^'!:]{#C)(xr„ - xom) = / dibld^k- (3.9) 

Jc 

Let us now, for the sake of simplicity, focus on the case of a defect line L which (i) 
can itself be represented as a single C", closed or not, parametric line without multiple 
points except possibly its extremities and (ii) is isolated in the sense that each of its 
points X is located inside a smooth surface S{x) bounded by a loop C{x) and such 
that S{x) \ {x} C Ql- Such a defect line L will be called an isolated dislocation or 
disclination. The jump [tol] of the rotation vector ajj around L is defined as the jump 
of ujI along #C(i) and hence is the same for any x and suitable C{x). Similarly, the 
jump [bl] of the vector bl around L is defined as the jump of b^ along #C(£) and is 
also the same for any x and suitable C{x), given xq. In fact, the following result is 
well-known (Kleinert, 1989): 

Theorem 3.1 [Weingarten] The rotation vector uj^ is a multifunction of index 1 on 
r^L whose jump Q% := [uJk] around L is an invariant of the defect-line L. Moreover, 
for a given xo, the vector b^. is a multifunction of index 1 on Ql whose jump B^ := [6^] 
around L is an invariant of the defect-line. 

Proposition 3.1 [Multiple-valued displacement field] ^From a symmetric smooth lin- 
ear strain tensor £*j on 0,l and a point xq where the displacement and rotation are 
given, a multivalued displacement field u* of index 2 can be constructed on Ql such that 
the symmetric part of the deformation gradient dj^^u* is the single-valued strain tensor 
£*j on Q,L while its skew-symmetric part is the multivalued tensor :— —eijkuj^. 

^We note that the curve C could be non rcctifiablc, i.e. of infinite length. In fact, integrals 
on fractal curves and the related Stokes' and Gauss-Green's theorems are analysed by Harrison 
and Norton (1992), where it is shown, by the C°°-smoothnoss of the differential form dmU)'!^dxm 
on that Eq. (3.7) still holds even when the Hausdorff dimension of C is higher than 1. 



8 



^From this result, the Frank and Burgers vectors can be defined as invariants of the 
single isolated line L. 

Definition 3.3 [Frank and Burgers vectors] The Frank vector of the line L is the 
invariant 

m ~ [loI], (3.10) 
while for a given reference pomt xq its Burgers vector is the invariant 

Bl ■- [hi] = [ul]{x) - eMm^*l{Xm - XQm). (3.11) 

A defect line with non-vanishing Frank vector is called a disclination while a defect 
line with non-vanishing Burgers vector is called a dislocation. 

Clearly a disclination should always be considered as a dislocation by appropriate 
choice of xo while the reverse statement is false since Ql might vanish. This is why 
in the present paper, the word "dislocation" means in the general sense a dislocation 
and/or a disclination. A pure dislocation is a dislocation with vanishing Frank vector. 

Remarli 3.1 Two distinct reference points xo and x'^ define two distinct Burgers vec- 
tors, related by 

Bk — B'i* = ekim{xom — x'o^)^ll , 

in such a way that B^Qk is an invariant independent of the arbitrary choice of xo- 
Therefore, for a non-zero Frank vector, the vanishing of the Burgers vector depends on 
the arbitrary choice of xq . 

The following result can be readily shown and is fundamental in the framework of our 
investigations since it implies conservation laws at the meso- and macro-scales. 

Theorem 3.2 Single disclination and dislocation lines are always closed or end at the 
boundary of Q,. Moreover, 

^^e*J = 

Definition 3.4 (Mesoscopic strain incompatibility) According to Eg. (3.1) com- 
bined with Eg. (3.3), the incompatibility tensor is defined by 

Vlk •= ei„-i„dmdnUJk- 

The strain field is called compatible on the set U if the associated incompatibility tensor 
vanishes on U. 

4 Distributional analysis of incompatibility for 
a single rectilinear dislocation 

4.1 The 2D model for rectilinear dislocations 

2D elasticity means that the strain £*j is independent of the "vertical" coordinate z. 
However this assumption introduces no restriction on the dependence of the multiple- 
valued displacement and rotation fields upon z. 
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Notations 4.1 Henceforth the single defect line will be assumed to be located along 
the z-axis. The two planar coordinates will be denoted by x and y or Xc,. The projec- 
tion of X = {xc, , z) on L is X — [0,0, z). By convention, Latin indices i, j, k,l, ■ ■ ■ take 
their values from 1 to 3 and are basically used for 3D elasticity, while Greek indices 
a, j3,^,5, - ■ ■ take the values 1 or 2 and are used for 2D elasticity. Symbols [e^, By, e^) 
or [ccCz) denote the Cartesian base vectors, while {er,ee,ez) denote the local cylin- 
drical base vectors. For a planar curve C, the notation dCa{x) = ea/sdxjs will be used 
for the curve normal. 

Let us observe that many fields are singular at the origin and tliat Ql is in fact 
tire domain where the laws of linear elasticity apply. Moreover, the strain can be 
decomposed into three tensors: 

^ij ~ SaiS/jjS^p -\- (^SizSj^S*, -It SjzSi^S*^'^ -\-5iz5jzS*z- 

planar strain 3D shear pure vertical compression 

The following propositions can be readily proved from Assumption 3.1: 

Proposition 4.1 [2D compatibility] In Ql, from 2D strain compatibility, there are 
real numbers K, and b such that 

!ea~,€0sdadp£*s = 0, 
eal3da£^z = K, (4.1) 
£^22 — 0(^Xci ~t~ b. 

Lemma 4.1 Let C{x) denote a family of 2D closed rectifiable curves. Then, in 2D 
elasticity, the Frank tensor and the strain verify the relation 

lim / Xadpujldx/s + e^pEp^dXa = 0, 

ic(4) 



C(i)-.i 



provided the length of C is uniformly bounded and as long as the convergence C{x) x 
is understood in the Hausdorff sense, i. e. in such a way that 

maxjlls — £11,1 G C{x)} — » 0. 
Proof. The second compatibility condition of Eq. (4.1) is equivalent to 

djSfiz ~ dfjS-iz = Ke^p, 
from which, so far as 2D elasticity is concerned, 

dpui^ := CK^d-fS^z = eK'ydf3£*z — KSk.p, 

and 

[xadi3LL>K + 5afteK~f£*iz) ~ dp (^Xaei.,j£*z^ — XcKSk^. 

Since, under the limit assumptions of this lemma. 



lim / XadxK — 0, 

and since the strain is a single-valued tensor, the proof is achieved. □ 
Lemma 4.2 In 2D elasticity the planar Frank vector Q.^ vanishes. 
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Proof. Since 



the planar Burgers vector simply writes as 

B* = [£^^ + eTj{x^ — xo-,)Spijj*) dxp - eT^{z ~ zo)iltf, 
Jc 

where C is any planar loop. By Weingarten's theorems the Burgers vector is a con- 
stant while the integrand is independent of z, from which the result obviously follows. □ 



In general, the present theory does not make any use of the linear elasticity constitu- 
tive laws and of the momentum and energy conservation laws, since in the framework 
of Continuum Mechanics arbitrary body forces and heat supply could be applied. 
However the remaining of this section will be devoted to present the three classical 
examples of 2D line-defects for which the medium is assumed to be body force free 
and isothermal (detailed computations are given in Appendix A). 

• Pure screw dislocation. From Eq. (A. 7), (A. 9), and Proposition A.l, the dis- 
placement and rotation vectors write as 



Uid = — and LOiCi = -V X Uid = —^Cr 

2tv 2 47rr 



(4.2) 



in such a way that the jump [uj*] vanishes identically, while from Eq. (A. 8) the 
Cartesian strain writes as 



-b: 



Moreover, in fii, appealing to Eq. 



y 







.3), the Frank tensor writes as 



47^7-2 



cos 26) 
sin 26 




sin 26 
- cos 261 




(4.3) 



(4.4) 



Pure edge dislocation. From Eq. (A. 6), (A. 7), and Proposition A.l, the dis- 
placement is the vector 

* -i?;(iogi + i) me 

while the rotation ui* vanishes together with its jump. The Cartesian strain 
writes from Eq. (A. 5) as 



[St 



-B*y 

27rr2 



V 







(4.5) 



noting that the tensor dmi^k = vanishes identically in Ql- 
• Wedge disclination. From Eqs. (A. 7), (A. 6), and Proposition A.l, the rotation 



vector IS 



with the multiple- valued planar displacement field given by 



(4.6) 
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and a vanishing Burgers vector: 

B* - iBl = [ul] - i[Uy\ + + ix) = 0. 
The Cartesian strain writes from Eq. (A. 5) as 



n^l - V*) 
in 



(log^ + l) 



and hence 



2nr 






cos 261 
sin 29 












(log ^ + 1) 




sin 29 
- cos 261 




sm 
— cos ( 




(4.7) 



(4.8) 



Remark 4.1 It should be noted that the energy density (or compliance) £* = ^a*j£*j 
is not L -integrable for both kinds of dislocations, while it is finite for the wedge discli- 
nation. Therefore, a Hadamard finite part distribution (Schwartz, 1957; Estrada and 
Kanwal, 1989) is needed to represent the compliance at the meso-scale (another ap- 
proach makes use of strain mollification by a so-called core tensor (Koslowski et al., 
2002)). This issue, whose solution requires to develop matched asymptotic expansions 
around the singular line in accordance with the infinitesimal displacement hypothesis, 
will not be addressed further in the present paper which only focuses on the geometry 
of dislocations. 



4.2 Mesoscopic incompatibility for a single defect line 

For 2D problems the incompatibility vector contains all the information provided 
by the general incompatibility tensor. The latter expresses on the one hand the non- 
commutative action of the defect line over the second derivatives of the rotation vector 
and on the other hand is related to concentrated effects of the Frank and Burgers 
vectors along the defect line. 

Definition 4.1 In the 2D case, the mesoscopic incompatibility vector is defined by 

nl ~ Ccpdadpujl. (4.9) 

A strain field is compatible if the associated incompatibility vector vanishes. 

As shown in the following sections, concentration effects will be represented by means 
of first- and second-order distributions. 

Notations 4.2 Recalling Notations 4-1, and stand for the sets SI, := {x £ 
fl such that X = (xa,z)} and Sl^ := f2z \ L, while the radius r = \\x — x\\ is the 
distance from a point x inside Q to L. Then, the ID Hausdorff measure concentrated on 
L is denoted by Sl (cf Ambrosio et al. (2000), Evans and Gariepy (1992) and Mattila 
(1995) for general definitions and properties on the geometric measure theory). 

In what follows the hypothesis consists in assuming that the strain radial dependence 
in the vicinity of L is less singular than a critical threshold. This is verified for instance 
by the wedge disclination whose strain radial behaviour is 0(ln r) and by the screw and 
edge dislocations whose strains are 0(r~^).^ For a straight defect-line L, according to 
these examples, the hypotheses on the strain and Frank tensors read as follows: 

■^A function /(e) is said to be O (g(e)) (e — ► 0+) if there exists K. tQ > s.t. < e < eo =^ 

l/(<:)| < K\g{()\- A function /(e) is said to be o{g{t)) (e 0+) if lim = 0. 

e^o+ 3(e) 
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Assumption 4.1 [2D strain for line defects] The strain tensor £*j is independent of 
the vertical coordinate z, is compatible on Ql in the sense that conditions (4-1) hold, 
is smooth on and -integrable on SI. 

Assumption 4.2 [Local behaviour] The strain tensor £*j is assumed to be o(r^^) 
(e 0^) while the Frank tensor is assumed to be o(r^'^)(e — > 0"*"). 

Theorem 4.1 [Main result for a single line] Under Assumption 4-1 and 4-2, for a 
dislocation located along the z-axis, incompatibility as defined by Eq. (4-9) is the 
vectorial first order distribution 

Vk = SkzVz + SkKvl, 

where its vertical component is 

ril = D,ISl + {B* - ei3jXoi3D,l) O^Sl, (4.10) 

while its planar components are 

= ^eKcBldaSL- (4.11) 

Proof. For some small enough e > and using Notations 3.1 a tube 0,, can be 
constructed around L and inside fl. Assuming that the smooth 3D test-function (f has 
its compact support containing part of L, Q^^^ denotes the slice of SI \ 0^ obtained for 
a given x £ L, i.e. 

fie.z {x G Qz such that ||a;Q|| > e}, 
while the boundary circle of Q^^z is designated by Ce,z. 

▲ Let us firstly treat the left-hand side of Eq. (4.10). Indeed, from Definition 4.1 with 
Eq. (3.2), Definition 3.1, and Eqs. (3.3) and (3.4), it follows that 



<r;fe,(/p>= / dz lim nk{z,tp,e) 



where 



Ilk{z,ip,e) ■— - / eal3dl3(^kda'pdS ~ / €apek~,„£pnda'fidCj. 

The boundedness of Idrdsipl on Q,l provides the following Taylor expansions of if and 
of daip around x: 

(p{x) = ip{x) + ruadaifiix) + Yi'Thdrdsip {x + ■yi{x - x)) , (4.12) 
daifiix) = daipix) + riyT-dT-ddip {x + j2ix - x)) , (4.13) 
with < 7i(a; — x), 72(2; — x) < 1. 

A Consider the first term of Eq. (4.12), noted Ilfc. By virtue of strain compatibility 
on SIl and Gauss-Green's theorem, this term writes as 



^k{z,ip,e) 



- / d~, {e^pOpLulf) dS = e~fi3dpujl<fidC~, 



Since by Notations 4.1 then Eq. (4.12) and Assumption 4.2 

show that, for e — > 0''', 

flk = J eji3di3Ujl(^(p{x) + Xada'p{x)'^dCj + o{l). 
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▲ Consider the second term of Eq. (4.12), noted 11^. On account of Assumption 4.2 
and by expansion (4.13), this term may be rewritten as 

Ul{z,ip,e) ■— — / eapek-fnSpndafdC-y 



▲ From Weingarten's theorem, recalhng that dC-y = e-yrdxr, the expression Ilfc = 
life + life writes as 

rife = (9a(/p(i) / (Xadri^k - ^aO^k-yn^lTEpn) dx^ 

+ Qlf{x)+o{l). (4.14) 

▲ Consider the first term of Eq. (4.14), noted life, and take 5 = 7 in the identity 

^kSn^-fT = Skz {SjsSriT — S„^6ts) — S„z {S^sSkr — Sk-iSrs) (4-15) 

in such a way that 



U'f.da<fi{x) / [XadTOjl - Skzfiap£pT + (^fer Eaflf ) dx^-. 



(4.16) 



▲ The cases k = z and k — k, are treated separately. 

• When k = z, Definition 3.2 shows that 

d/sb* ~ £p^ + eTj{xj - xo^)d0Ujl — erj{z — zo)di3Lo* 

which, after multiplication by e^a and using Eq. (4.15) with r, a and z substi- 
tuted for k,S and n, is inserted into Eq. (4.16), thence yielding: 

= dc,tp{x) / {eradfjbl + XQadfjUjl + (2 - zo)di3ijjl) dxfi, (4.17) 

and consequently, from the definitions of the Frank and Burgers vectors, 

lim n'^ = < {ecT-B* - {z - zo)Q.a - ajoa^t} 9a(5o, ifz », (4-18) 

where So is the 2D Dirac measure located at and ipz{xa) := (p{xa,z), while 
symbol <C •, • 3> denotes the 2D distribution by test-function product. 

• When k — n, Definition 3.2 shows that 

dpbl ~ + ejT{xj - xoj)dpLO*, 
from which, after multiplication by e^a, it results that: 

Xadrl^K = -e^adrbl + €Kcy£^z + XQadrLO* + {x^ - XOK)dr(^i- 

Then, by Lemma 4.1 with a permutation of indices k and a, Eq. (4.16) also 
writes as 

n'„ = da<p{x) I {-e^adisbl + e^aS^z + XOadfliiJK — Xo^dfjUJa) dx,3 

-fo(l). 
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On the other hand, from Eq. (4.16) and Lemma 4.1 (i.e. from strain compati- 
bility) it follows that: 



n' 



^a^pix) / [-e^pE^^dXa + eapSpzdxK.) + 0(1) 

^afix) eaK£^:,dX0 + 0(1). (4.19) 



By summing this latter expression of n'„ with Eq. (4.19), from the definitions 
of the Frank and Burgers vector it follows that 



n'„ = -da(p{x)eaK {B* — €j/3Q*xop) + 0(1). 
Hence, in the limit e — > 0"*" Eq. (4.20) writes as 

lim II'k =< I ^e^aB* - ^e^ai^^fjQtjXofi I daSo, (pz > 

Therefore, the result is proved on since 

lim Ilk{z,(fi,e)= lim Il',,{z,tp,€)+ ^ illSo, (fiz ■ 



(4.20) 



(4.21) 



(4.22) 



▲ As suggested by Eq. (4.12), to obtain the result for the entire domain Q it suffices 
to integrate Eqs. (4.17) and (4.20) and expression illip{x) over L, in order to replace 
So by the line measure 5l in Eqs. (4.18), (4.21) and (4.22). By Eqs. (4.12), (4.18), 
(4.21) and (4.22), the proof is achieved. □ 



4.3 Applications of the main result 

Throughout this section {x,y,z) denotes a generic point of f2_L and all tensors are 
written in matrix form in the Cartesian base {ex,ey,ez). 

• Screw disclocation. Since B* = fi* = 0, Eq. (4.11) yields 



dySL 





This result is easily verified with use of Eq. (3.6). One needs to compute 

/ ^kpn^a/sSfindpda^dV , that is to calculate the integral of 
Jn 



Bl 



By integration by parts, using Gauss-Green's theorem on fi, and recalling that 
test-functions have compact supports and that 9m log r = — ^, these integrals 
become 



Blf 



-9.^(c>.^+a,2i^) 





dV^ 



-dy(fidf-,Jogr 




dV. 



Hence, from A (log r) — 2it5l , the first statement is verified. 
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• Edge dislocation. Whereas dm^^k identically vanishes on fii, it is easily seen 
that Eqs. (4.10) and (4.11) with B* = fi* = yield 








We must compute [ryl.] = / epni^ea/iSpn^pdaipdV . For k = 1 and 2 and with 
Jn 

n 7^ 3, the tensor eap£fj„dpda'P equals Ey^dzdyip — Eyydzdxf and £*ydzdx^ ~ 
Sxxdzdyifi respectively. By integration by parts, the related integrals vanish. For 
fc = 3, the integrand is 

epnzeal3£p„dpda(p£xxdydyip + Eyydxdxif - 2£xydydx(p. 

Integration by parts provides the expression / -dx(pA{logr)dV, achieving 

in 27r 

the second verification. 



Wedge disclinaUon. Incompatibility reads 

IVk] = ^ 






5l 



We must calculate < r/l, (p >. For k — 1, k = 2, n^O and p — 3 the integrand 
vanishes. For fc = 3, one computes 

o ^ f^z(l — I'*) A /n r. fl*z(l — v*) . r. 

epnei^Smndpdiip = '-fA{\og—)A — '-If A{\og—) 

Aiv H in H 



4n 



achieving the third verification. 



5 Distributional analysis of incompatibility for 
a set of isolated dislocations 

In the previous section, a single defect line was considered. However, to address 
the macro-scale physics, homogenisation must be performed on a set of dislocation 
lines whose number tends to infinity in order to define regular defect density tensors. 
Therefore, our goal is to introduce appropriate hypotheses that can easily be applied 
to a set of defect lines and to a regular defect density as well. 

5.1 Governing assumptions for the strain and Prank ten- 
sors 

Besides the strain Assumptions 4.1 two measure hypotheses on the strain derivatives 
are introduced to replace the local Assumption 4.2 in order to validate Theorem 4.1 
in a global framework. 

Assumption 5.1 The strain divergence and trace gradient daSai and d-fS^K '''Te finite 
Radon measures on Q.. 

*A finite Radon measure on Q is a measure bounded on compact subsets of C 
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Remark 5.1 No assumption could be made on the complete Lebesgue integrable strain 
gradient without contradicting the 2D examples of Appendix A. On the other hand, it 
can be shown that the sharp Assumption 5.1 are required to demonstrate Proposition 
5.1 

Remark 5.2 Assumption 5.1 is natural in infinitesimal elasticity if one considers the 
strain-stress constitutive law and the equilibrium laws. As a consequence, the stress 
divergence must be a measure on O,. 

The following Lemmas are needed for the proof of Proposition 5.4. 

Lemma 5.1 • A solenoidal distributional vector field on Slz writes as 

Ua = £a-/d^(j), (5.1) 

with (j> eV'{fl^). 

• A symmetric solenoidal distribution tensor a^p on Q,^ writes as 

a-afi = eaj£f}TdjdTi>, (5.2) 

with ip ev'in^). 

Proof. 

• First statement. Let (t>o be any 3;2-primitive distribution of ai (Schwartz, 1957). 
Then 92 0o = cti and, from the solenoidal property of Oq,, there exists a distri- 
bution G{xi) s.t. ditjjQ + 0,2 = G{xi). By a::i-primitivation of G{xi), it is easy 
to find F{xi) s.t. diF — G{xi), and to verify that (f) = (j)o + F{xi) solves the 
problem. 



• Second statement. /,From daaap = 0, there is a distribution (j>i3 s.t. aap = 
tajd-f4)p. Then ea-jd-, {dpcjiiB) dpaap = and hence dpcfip is a constant G or 
equivalently 9/3(0/3 — ^Gxp) = 0. From Eq. (5.1), there exists a distribution ip 
such that 0/3 — ^Gx0 — eprdr'ipj a-nd hence Oap = ea-fCj^rd-^dTip + |ea/3C. The 
symmetry of aafi implies that G = Q. □ 

Lemma 5.2 • For a given V' {Q, 2) -scalar function f, there exists an irrotational 
distribution field such that 

dpgp = /. (5.3) 

• For a given (flz)-vector function fp such that dpf/s — Ag where g is a (fiz) 
function, there exists a symmetric compatible tensor gap on such that 



dagap = //3. (5.4) 



Proof. 



First statement. It is sufficient to consider an ultra-weak solution (Brezis, 1983) 
of AH — f and to define gp = dpH. 



• Second statement. By primitivation, there is a non-compatible {^lz)-field g^p 
such that /i = 9iPii,/2 = 92<?22 a-nd = (?2i ~ 5i2- A necessary condition 
for gap to exist is that gap = gap — gap verifies dafjap ~ 0, or by Lemma 5.1 
that gap = ea-yepTd-,dT4> for some gauge distribution cf). In order that gap be 
compatible on flz, (f> must satisfy the following equation, equivalent to the 2D 
compatibility of gap on fi^: 

AA0 = AffL - dpdaQap = A(5L - g)- (5.5) 
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Up to a harmonic and hence smooth function on Q^, the solution of Eq. (5.5) is 
the solution of A(/) — p*^ — g. Since the right-hand side is L^{^lz), a solution 
exists in the ultra- weak sense and hence the existence of a symmetric compatible 
distribution field g^p on fiz verifying Eq. (5.4) follows. □ 

Lemma 5.3 For constant C and Ci3, there are a vector g^ and a symmetric, compat- 
ible tensor Gctfi on such that 

d^g^ = CSo, (5.6) 

daGaP = CfjSo. (5.7) 

Proof. The solutions are given by g^ = (27r)~^9„ logr and Gap ~ ^ {dcHp +dpHa), 
where 

= ^(hoer-^]-C,^^''' 



□ 

Lemma 5.4 Under Assumptions 4-1 and 5.1, the strain components can be put in the 
form; 

£*, = E^ + e^, (5.8) 

£ap ~ Eap + Ca/J, (5.9) 

where vector has a vanishing curl on for any given z while vector is o(r^^) 
as r —>■ 0^ , and where tensor Eafs is compatible on for any given z while tensor 
ecB IS o{r~^) as r — > 0'''. 

Proof. By Assumption 5.1, OkSki is a Radon measure on ^1., and hence writes by 
Radon-Nykodym's decomposition theorem as 

d^s:, = f, + (l>,, (5.10) 

where fi £ L^{fl~) and where (j>i is a Radon measure on Q~ singular with respect to 
Lebesgue's measure. As a mere consequence of the smoothness of OkEki on 4>i is 
a concentrated measure on 0.^ and hence is proportional to the Dirac mass 5o, 

0> = C,5o = (27r)-'aa2logr. (5.11) 

▲ First statement. 

• By Eqs. (5.10), (5.11) with i = z, and Lemma 5.2, there exists an irrotational 
g„ such that 

{£L - 9. - {27v)-'C,d^ logr) = 0, 

in such a way that, by Lemma 5.1, 

- ,g« - (27r)~'a9«logr = e«^9^^, (5.12) 

where tp is a, distribution. Apply the curl operator to Eq. (5.12) and take into 
account the irrotational property of g^ in such a way that Ai/) = ei^^di3£^z. 
Since S^z is a L^-vector, its curl is a first-order distribution^ and hence, by the 
strain compatibility which ensures the curl of E^z to be a constant K on $1° 
and a combination of the Dirac mass and its first-order derivatives at the origin 
(Schwartz, 1957), writes as K + cS + c-yd~,5. 



^Following Schwartz (1957), a distribution is of order 1 if it defines a Hnear continuous map 
onCi(n). 
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• Now in the resulting equation 

dp (^e^pS^z - (27r)"^ca^logr- - y^-fi^ = Cjd^S, (5.13) 

the term on the left-hand side is the divergence of a L^-vector, and hence Eq. 
(5.13) has no distributional solution unless c-y = 0. 

• It results that Atp — K + (27r)~^cA(log r) provides a gauge field ip which writes 
as 

= /i + (27r)"^clogr, (5.14) 

where /i is a solution of Ah = K on Q^. It is easily verified that the curl of tp is 
o(r~^) as r ^ 0"*". 

• Defining = + (27r)~^Cz9K log r and = ^K-ydjip in Eq. (5.12) achieves 
the first statement proof. 

▲ Second statement. 

• Let us prove that the divergence of ft is the Laplacian of an L^{fl:,) function. 
In fact, since rj'^ writes as 

rj* = da {da£^^ - di3£*p) , (5.15) 

it is from Assumption 5.1 a concentrated first-order distribution writing as a 
combination of the Dirac mass and its first-order derivatives. Hence: 

9/3 ff3 = dadpEap - dpcj>p = A£*^ -V*z- dp(f>p = - c5o - Cjd^,So 

= A (f*« - clogr - c^9^1ogr) , (5.16) 

where c,Cj,c,c-f are constants. 

• From Eqs. (5.10), (5.7), (5.4) and Lemma 5.2, there exists a compatible g^p 
such that 

d^{£:p-g,p-G^p) ^0, (5.17) 

in such a way that, by Lemma 5.1, 

£kP - 9kP - Gk/3 = e^^ep^djdrA, (5.18) 

for some gauge field A G D'{Q,z) verifying, by the compatibility of g^p and G^p 
on Q,z, the relation 

rj* = AAA on Q,. (5.19) 

Hence, since the left-hand side writes as a combination of derivatives of So of 
order lower or equal to 1, the field A is the solution of AA = {a + a-yd-y) logr 
with constant a, a-,, up to a smooth harmonic function on Q.^. It follows that 

A — {a + ajdj) ^^(logr — 1)^ is a C"{Qz) solution of Eq. (5.19) such that: 

d^dpA is o(r"^) as r ^ 0+. (5.20) 

• The proof is complete with the definitions Ek,p — Gnp+g^p and e^p = ei^^ep-rd^^idrA 
in Eqs. (5.18) and (5.20). □ 
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5.2 Mesoscopic incompatibility for a set of isolated defect 
lines 

Theorem 5.1 [Main 2D result] Under Assumptions 4-1 o,nd 5.1, for a set C of iso- 
lated dislocations parallel to the z-axis and located at the positions Xp, L G C, incom- 
patibility as defined by Eq. (4-9) is the vectorial first order distribution 

Vk = Skzr]l + Sk^nl, (5.21) 

where 

• its vertical component is 

77* = ^ (Q.tSL + £07 (^B* + e/3^,{xp - xop)Ol^ daSh^ , (5.22) 
Lec 

• its planar components are 

77k = ^ ^eKaB*da5L- (5.23) 

Proof. ^From Lemma 5.4 the strain (n = a or z ) is decomposed in compatible 
parts {E/s and Eap) and o(r~^) parts (e/3 and Bcb) to which the demonstration may be 
limited by linearity. Since from Eqs. (5.14) and (5.20) the gradients d-yCp, d-yeafi axe 
o{r~^) for r — > 0^, the proof of Theorem 4.1 can be followed for every L £ C as soon 
as is replaced by ep and f^^,. by ep-r- However, since the dislocations are located 
at positions a;^ instead of 0, an additional shift 2:^ is required in Eq. (5.22). □ 

5.3 Mesoscopic defect densities in 2D incompatible elas- 
ticity 

Since the tensors 0*^., A*j,, vanish for i^ z, the 2D densities for an ensemble C of 
rectilinear dislocations write as follows'': 

Definition 5.1 

e: — ^4.nt^<5L (5.24) 

At. ~ ^B:^<5l, (5.25) 

Lec 

Qfc alk = Al — Skaeai3Q*,{x0 — X013). (5.26) 
Moreover, in the 2D case, the contortion tensor writes as: 

= Siza* — -a*,5ij. (5.27) 
The following result expresses the incompatibility in terms of : 

Theorem 5.2 Under Assumptions 4-1 and 5.1, the mesoscopic strain incompatibility 
for a set C of rectilinear dislocations writes as 

Vk = ©fc + eal3daK,lf3, (5.28) 

or equivalently as r/^ = + ekaidaK,*i. □ 



^Various notations arc used in the literature to represent the defect densities. In particular, 
Nye (1953), Kroner (1980) and Kleinert (1989) give different definitions of the dislocation 
density and contortion tensors (without considering disclinations in the first two cases). We 
here follow Kroner's and Kleinert's notations for q*^ and Nye's original definition of k*^, with 
Nye's a* J here denoted by a*^ . It should be recalled that the term " contortion" was introduced 
by Kondo (1952). 
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Proof. Consider any straight dislocation L ^ C located at a given G Q.. From 
Theorem 5.1, incompatibility writes as 

nl = Skz {^15l + ea-i [B* + €is^{Xis - Xois)Q,*,) da5L) 

+ 5kn^e^cBldc,5L. (5.29) 
Taking into account Eqs. (5.24), (5.25), (5.26), and (5.27), and the relation 

da {{Xp - Xoi3)dL) = da (^{Xp — a;o/3)5i^ = {X0 - Xop)da5L, 

it results from Theorem 5.1 that incompatibility can be written in the alternative 
formulation 

■nl{x^) = Ol{x^) + eapdai^lpix^), (5.30) 

or equivalently as vj^ix^) = Q'^ix^) + ekaidai^ziix'^)- The result follows after summa- 
tion on L e £ and using Eqs. (5.24), (5.25), (5.26), and (5.27). □ 

First of all, the tensor djdml is defined on the entire f2 in a similar way as djujl: 
Definition 5.2 

djdiul ~ djSki + f^kpidjujp. (5.31) 

By Proposition 3.1, the displacement field u^. is a multivalued function of index 2, 
which is obtained on Ql by recursive line integration of Sj^'Sj^^'u^ — dj"'' {£^1 +'^fci) 
and hence by recursive integration of djdiul. 

Remark 5.3 In the situation where, for a particular selection of the reference point, 
the dislocations have vanishing Burgers vectors, the disclination density equals the 
incompatibility 

tapdadptol = Ql = ril. (5.32) 

Using an arbitrary reference point, this expression is certainly false in the general case 
where disclinations coexist with dislocations. Moreover, the tensor djdiu^ does not 
provide relevant information in terms of defect densities since eijidjOiu^ = on 

The mesoscopic vectors and tensors 0^,A^,a^,K^ and r]l are concentrated distribu- 
tions on the defect lines which provide all the information on dislocation and discli- 
nation densities. However, homogenisation to the macro-scale still requires to clarify 
their link with the multiple-valued rotation and displacement fields. In order to re- 
solve this problem, the tensors djUjl and djdiul are completed by appropriate con- 
centrated efi^ects in the defect lines, without however modifying their relationship with 
the multiple- valued displacement and rotation fields defined in S^l. 

Definition 5.3 

Bpujl := dfjujl - (5.33) 

9a3/3Wfe := dadl3ul — ekp0lipa = da£kl3 + <^kp0<iai^p- (5.34) 

Theorem 5.3 The vector and tensor distributions S/joj^ and QadpuX verify: 

MESOSCOPIC DISCLINATION DENSITY Ql = €apda5i3iol, (5.35) 

MESOSCOPIC DISLOCATION DENSITY al ^ ea/sSaQpul- (5.36) 
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Proof. The first statement is a mere consequence of Eq. (5.28) while the second one 
follows from Eq. (5.33) by simple calculations, noting that dmi^m = on i7 and that 

"fc = i^lk - i^lpSzk- (5.37) 

□ 

Remark 5.4 Eqs. (5.35) and (5.24) directly show that 

I ec^pd^dpuoldS = I eldS= J2 "fc"^' (5.38) 

where the domain S is bounded by the counterclockwise-oriented Jordan curve C, which 
encloses once each defect line of the subset Cc of C. Similarly, Eqs. (5.36) and (5.26) 
show that 

/ eaffScQpuldS = / {Al - 5ka<^a,BQ*z{x0 - Xoi3))dS, 

J s J s 

= 5Z (^Bl^ - Skae^ii^l^ {x^ - xoisij . (5.39) 

Remark 5.5 The vector diuj* does not verify Stokes theorem, neither in the classical 
sense, since eaftdad pijj*^ is singular at , nor in a measure theoretical sense, since 
eai3dadi3Lo* IS not a measure but a first-order distribution given by Eq. (5.22). As 
often observed in the literature, even in an inappropriate context, a formal use of 
Stokes theorem may give a correct final result. We here prefer to avoid any confusion 
and hence to mention that, in view of a clarification of Stokes ' theorem in the context 
of defective crystals, the following formula holds and can be proved as a consequence 
of the previous definitions: 

i dlLuldXl = / €al3dalilijJ*zdS. (5.40) 
Jc J Sc 

6 Macroscopic analysis 

6.1 A first approach to homogenisation from meso- to 
macro-scale 

The mesoscopic results given in the previous sections are now homogenised (in an ap- 
propriate manner, whose description (Kroner, 2001) is not the purpose of this paper). 
Indeed, in the context of linear elasticity, the macroscopic elastic strain £ij is obtained 
by averaging the mesoscopic stress a*j and hence the macroscopic elastic incompati- 
bility rjik is obtained by averaging the mesoscopic incompatibility ry*;.. Moreover the 
defect densities are homogenised and the macroscopic counterparts of 0^, AJJ;,afe and 
write as Ofe,Afe,Qfc, and Hij, with 

ctk = K,zk — i^ppSzk and Kij — Si^aj — -a^Sij . (6.1) 



Definition 6.1 (Macroscopic Frank and Burgers tensors) The Frank and 
ers vectors crossing a macroscopic surface S are defined as 



i^kiS) := / OkdS, (6.2) 
Js 

BkiS) := f AkdS. (6.3) 
Js 
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By homogenisation of Eqs. (5.33) and (5.34), the macroscopic counterparts of Defini- 
tion 5.3 write as follows: 



Definition 6.2 

dpUJk ■■= tkpqdpEqP - Kfc/3, (6.4) 
SaS/S^ifc := daSkl3 + <^kp0(^atOp, (6.5) 

where £kp and Kkp define the macroscopic elastic strain and contortion. 

Moreover, the macroscopic counterpart of Theorem 5.2 (i.e. the fundamental equation 
"inc £ = Q + curl k" of the continuum theory of defects by Kroner (1980)^) and 
Theorem 5.3 now follow from homogenisation of the mesoscopic defect densities and 
from Definition 6.2: 

Theorem 6.1 

KRONER'S IDENTITY nk = Qk + eafidaHkil-, (6.6) 

MACROSCOPIC DISCLINATION DENSITY Qk = eapda(i fii^k, (6.7) 

MACROSCOPIC DISLOCATION DENSITY ak ^ e^fj^a^Q fiUk- (6.8) 

Remark 6.1 By Stokes' theorem, if S is a region enclosed by a curve C , which might 

have only fractal regularity (Harrison and Norton, 1992), then Q,k{S) — I dfsLUkdx/}. 

Jc 

Moreover, m the absence of disclinations, Bk{S) = / OkdS. 

Js 

The macroscopic density tensors and ftij, as obtained from the single- valued meso- 
scopic densities, have a geometrical interpretation (Kroner, 1980; Anthony, 1970) 
which will be discussed in the following section. Indeed, ak is directly related to the 
torsion of a body submitted to an incompatible purely elastic deformation to which a 
non-Riemannian connexion is attached due to the contortion Kij. 

6.2 The non-Riemannian macroscopic body 

The following geometric objects are introduced after homogenisation of the well- 
defined mesoscopic elastic strain and defect densities, in order to provide the model of 
a macroscopic body endowed with a law of parallel displacement together with internal 
torsion accounting for the defective crystal structure. 

Definition 6.3 

METRIC TENSOR: Qij := 5ij - 2£ij, (6.9) 

TORSION: Tk-ij := —^<^pij (ctpk - ekmnQpm{Xn — Xo,i)) , (6.10) 

SYMMETRIC CHRISTOFFEL SYMBOLS: 

f fe;ij ■- i [diQkj + djQki - dkgij) , (6.11) 

CONTORTION: /\Vk■,^j := T^jk + Tiijfe - Tk.ji, (6.12) 

NON SYMMETRIC CHRISTOFFEL SYMBOLS: 

^k:ij ~ ffe;ij - Arfe;ij. (6.13) 

Remark 6.2 The metric of the actual configuration R{t) isSij. Therefore, as required 
( cf Introduction and Remark 2.1) the reference configuration Ro is nowhere used to 
introduce the above objects. 



'^Note that different sign conventions for the rotation vector and incompatibility apply in 
Kroner's work. 
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Since small displacements are considered, no distinction is to be made between upper 
and lower indices. 



Lemma 6.1 The tensor Qij defines a Riemannian metric. The symmetric Christoffel 
symbols Fk^Aj define a symmetric connexion compatible with this metric, while Tk-ij and 
Arfe;ij are skew-symmetric tensors w.r.t. i and j and i andk, respectively. Moreover, 
the components of Tk;ij for i — z or j = z vanish in the 2D case. 

Proof. The first statements follow from basic definitions (Dubrovin et al., 1992; 
Schouten, 1954) while the last one follows from the fact that, in the 2D case, Qpfe(a:^) 
and Opm{x'^) for L G C are proportional to Tp5l{x^) with Tp standing for the tangent 
vector to the defect line. □ 

Proposition 6.1 The Cristoffel symbols Tk-ij define a non-symmetric connexion com- 
patible with the metric gtj and whose torsion writes as Tk-ij. ^ 

Proof. It is easy to verify (Dubrovin et al., 1992) that Tk-,ij is a connexion since 
Tk-,ij is a connexion and Arfe;^^ is a tensor. Denoting by Vfc (resp. V^) the covariant 
gradient w.r.t. Vk-,ij (resp. Tk-,ij), and recalling that a connexion is compatible with 
the metric Qij if the covariant gradient of gij w.r.t. this connexion vanishes, we find 
by Eq. (6.13) 

'^kQij '■ ~ dkQij — Tl-ikQlj — ri-jkSli 

= Vfc5y + AVi-ikgij + Ari-^,k9ii, (6.14) 

where in the right-hand side, the 1"* term vanishes by Lemma 6.1 while the 2"'' and 
3'^'' terms cancel each other since Ari-,jkgii ~ ATi-jk = —AFj-jk- It results that the 
connexion torsion, i.e. the skew-symmetric part of ATj-^k w.r.t. i and k, writes as 

i (Ar,;,fe - Ar,;fe,) = -i (Ar,;,fc - A^k-,,^) = i ((Ar^;,, - Ar,;,,) + 

(AFfc;,, - ATk-aj) - (Ar.;,fe - Ar,;fe,)). (6.15) 

Observing that the 1"' term in the right-hand side of Eq. (6.15) writes as Ark-,ij while, 
by Definition 6.3 (Eq. (6.12)), the left-hand side and the two remaining terms of the 
right-hand side of Eq. (6.15) are equal to Tj-ik,Tk-ji and —Ti-jt, respectively, the proof 
is complete. □ 

The following result shows AF k-,ij as directly linked to the contortion Hij. 

Proposition 6.2 In the 2D case, the contortion tensor Ark-,ij writes in terms of Kij 
as 

ATk-^ij — 5kK {SiaSjpeKalizfi) + SiaSjz£aTK.TK + &iz& j pepr I^tk 

— 5l^z&ia5jpeafiK,zz. 

Proof. For fc = 2, by Definition 6.3, the last statement of Lemma 6.1, and Eq. (6.1), 
it is found that AT^-^ij — AT^-apSiaSjis, with 

AFzjc/J — T^-^aP = ——CapO-z = —Capl^zz 

1 c 

— —-eaTOTfjOtz = ^arl^Tp. 



° In the literature, a so-called Bravais' crystal is a macroscopic body endowed with a lattice 
where parallel displacement along the crystallographic lines is defined by the connexion Tf^-ij 
of Theorem 6.1 and where the metric is not defined by Eq. (6.9), but by the motion of an 
internal observer who would measure his own displacement by counting the atomic lattice 
steps, without feeling the body torsion (Kroner, 1980). 
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For k — K, hy Definition 6.3 and the last statement of Lemma 6.1, it is found that 

with Tz;^K ~ e^T^TK and T^;ti, ~ ~^(.tv + e^.^Qz{x^ — 2:07)). Since the combina- 
tion of the terms in O, vanish in AVK\ii; the proof is completed by observing that 

The following definition introduces two differential forms related, on the one hand (by 
Definitions 6.2 and 6.3, and Theorem 6.1 and Proposition 6.1) to the homogenisation 
of the well-defined mesoscopic defect measures and, on the other hand, as shown by 
the forthcoming theorem, to macroscopic incompatible rotation and distortion vectors. 



Definition 6.4 The following differential forms are introduced: 

dujj := dpLjjdx/s, (6.16) 
dPu ■- -Ti,updxp. (6.17) 

In the literature the existence of an elastic macroscopic distortion field is generally 
postulated (Mura, 1987; Head et al, 1993; CermeUi and Gurtin, 2001, 2002; Koslowski 
et al., 2002; Ariza and Ortiz, 2005) and the global distortion decomposition in elastic 
and plastic parts follows'"' . The point of view of the present paper is to avoid this kind of 
a-priori decomposition, which we believe cannot be rigorously justified. Nevertheless, 
the following theorem introduces rotation and distortion fields (which are not the 
global rotation and distortion related to the macroscopic strain) in the absence of 
disclinations. In contrast with the classical literature where it is basically postulated 
that dislocation density is the distortion curl, this relationship is here well-proved. 

Theorem 6.2 [Bravais rotation and distortion fields] If the macroscopic disclination 
density vanishes, there exists rotation and distortion fields defined as 

rx 

BRAVAIS ROTATION u!j{x) := L0° + j duj, (6.18) 

J xo 

rx 

BRAVAIS DISTORTION I3u{x) := £m{x°) ~ eMjLO° + j dfiu, (6.19) 

J xo 

With pki ~ £ki — ^kijtj^j, and where uj*] is arbitrary and the integration is made on any 
line with endpoints xo and x. Moreover, 

daPkfi = SaS^Ufe and eapdaPkp = Ofe. (6.20) 

Proof. By Definition 6.3, the symmetric part of the connexion writes as 

— ^(i;k)i3dxp = -^dpgkidxp = -^dmgkidxm = dm£kldXm = d£kl, 
while, by Definition 6.3 and Proposition 6.2, the skew-symmetric part writes as 

— ^\i;k]f> ^ ~^{dkgii3 ~ digkfi) + /^Ti:k0 ^ dk£ii3 - di£ki3 + Ari-ki3- 

Observing, by Definitions 6.2 and 6.4 and Proposition 6.2, that dujj = dpujjdx/s 
= —^^ikj^[i;k]i3dxi3, it results that dPki = d£ki — ekijdujj. Under the assumption 
of a vanishing macroscopic disclination density, the existence of well-defined Bravais 
rotation and distortion fields follows from Eqs. (6.16) and (6.20), Theorem 6.1, and 
Remark 6.1. Moreover, since daPkp ~ da£kp — ekiij^Qai-Oj , by Eq. (6.5), it equals 
9ct9/3Wfc, completing the proof by Eq. (6.8). □ 



^In fact, the distortion is often considered as a constitutive variable in dislocation models 
(Davini, 1986; Gurtin, 2002; Ariza and Ortiz, 2005). 
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Remark 6.3 Referring to "Bravais" instead of "elastic" rotation and distortion fields 
is devoted to highlight that these quantities do not have a purely elastic meaning 

Remark 6.4 The Bravais distortion does not derive from a Bravais displacement in 
the presence of dislocations. In fact, around a closed loop C , even if the macroscopic 
disclination density vanishes, the displacement differential as defined by duk :— Pkadxa 
verifies by Theorem 6.2 the relationship: 

[ duk= [ ep^dpPkcdS = ak[S). (6.21) 
J c J s 

Remark 6.5 Eq. (6.18) indicates that symbol 9 m Eq. (6.16) becomes a true deriva- 
tion operator in the absence of disclinations. 

Remark 6.6 Theorem 6.1 defines an operation of parallel displacement according to 
the Bravais lattice geometry. The parallel displacement of any vector Vi along a curve 
of tangent vector dx^^ is such that dXa'^VaVi = and hence that the components 
of Vi vary according to the law d^^'^Vi = —^i-jpVj dx^l^^ (Dub rovin et ai, 1992). This 
shows the macroscopic Burgers vector and dislocation density together with the Bravais 
rotation and distortion fields as reminiscences of the defective crystal properties at the 
nanoscale. In fact, if dxi}\ dx^^-' are two infinitesimal vectors with the associated area 

dS := ev^dx^}'' dx^^\ it results from Eq. (6.10), Remark 6.1, and the skew symmetry 
ofTk-ai3 that, in the absence of disclinations, 

dBk = akdS = —Cap^kiapdS = -Tkia/sidx^;^ dx'jP - dx'jj^ dx'-^^ ) , 

whose right-hand side appears as a commutator verifying the relation 

dBk = Cap^a^fjUkdS = -eai3d^"\dx^'^^). 

7 Concluding remarks 

In this paper we have developed a 2D theory to analyse dislocated single crystals at 
the meso-scale by combining distributions with multiple-valued kinematic fields. The 
distributions are basically concentrated along the defect lines, which in turn form the 
branching lines of the multivalued fields. As a consequence of this analysis, a basic 
theorem relating the incompatibility tensor (as derived from the deformation field) to 
the Frank and Burgers vectors of the defect line has been established. This theory 
provides a framework for the homogenisation of the medium properties from meso- 
to macro-scale. In particular the macroscopic dislocation density is defined without 
stipulating an a-priori distorsion decomposition into elastic an plastic parts (which 
does not exist, actually). The classical relationship between Bravais distortion and 
dislocation densities, instead of being a definition, now appears as a result taking its 
origin from the meso-scale analysis. Moreover, the torsion and contortion tensors, 
which both describe the defective macroscopic crystal, are now properly understood 
as averages of concentrated mesoscopic tensors. Since the latter are the differentials 
(in an appropriate sense) of multivalued mesoscopic fields, we have shown how meso- 
scopic multivaluedness is recovered in the geometric properties of the non-Riemannian 
macroscopic crystal. In particular, in contrast with the mesoscale (where defects are 
due to the multivaluedness of the rotation and displacement fields) the macroscopic 
Burgers vector now appears as the commutator of a non-closed differential operator 
related to the body torsion. 

Extension to the 3-D case is under investigation. Here, the handling of non-rectilinear 
curves will be required in the framework of the geometric-measure theory. This should 
eventually make it possible to consider a set of defect curves, freely occupying the crys- 
tal with possible intersections and accumulation regions-forming so-called dislocation 
clusters. 
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A Appendix: Computation of 2D rectilinear dis- 
locations 



A.l First group of solutions: planar displacement field 

^From the constitutive law a*j = XE^^Sij + 2/ifj* with A,/i the Lame coefHcients and 
since £1^. — f*^, the following planar law holds; 

0"a/3 = hi*£*^SaP + 2/if (A.l) 



with the planar compressibility modulus k* defined by k* := A + ^ and the planar 
dpo*^^ = it follows that 



deviatoric strain given by f — £*0 — ^£*^Sai3- ^From the equilibrium conditions 



o"q,3 = ea^eesdjdsF, (A. 2) 
for a smooth enough Airy function F, in such way that 

a*^^^dlF = AF. (A.3) 
The relations between stress and strain are 



= (^ll3 - -^^j-ySalS, (A.4) 



with the 3-D and planar elastic coefficients E — ^^^^ — — — —, v — —- r, E* := 

A + ^ ' 2(A + ^)' 

E J » u 
-, and V ■- . 

1 — 1/^ I — u 

The first compatibility condition Eq. (4.f) writes from Eqs. (A.3) and (A.4) as 

AAF = 0. 

In this and the following sections, functions of the complex variable Z = x -\- iy and 
its conjugate Z are now introduced. Remembering that, compared to holomorphic 
functions, analytical functions may be multivalued it is easily seen that given two 
analytic functions / and g, all real functions of the form 

F = ^{Zf + g} 

satisfy Eq. (A. 2) and vice-versa. Eq. (A.3) then shows that 

ol, + oly^AW:{f{Z)}. 

/.From Eqs. (A.3) and (A.4) the deformation tensor is given by 



= ^^K{/'(Z)}, 



(A.5) 



f - + 2i£ly = 2Mp(Z/"(Z) + g"[Z)), 

yielding after integration 

i5*«-z<) = {i-v*)f{Z)-{l + u*)(Zf{Z)+g'[Z)), (A.6) 
E*u: = 43{/'(Z)}. 

It should be emphasised that f must be a single- valued field. 
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A. 2 Second group of solutions: vertical displacement field 

Another solution concerns the particular case where u* = 0, noting that each solution 
of 2D elasticity can be decomposed into a purely planar and a purely vertical solution. 
In fact, since stress equilibrium shows that 

(A + fi)didjU* + fiAu* = 0, 

it is easy to infer for i = z that 

= ^^^mnz)}, (a.7) 

with h{Z) an analytic function. Then 

2E 

The function h'{Z) must be uniform. The complex rotation is 



£L - iSy. = ^^^h'{Z). (A.8) 



2^ h'{Z). (A.9) 

In 2D isothermal linear elasticity without body forces, every displacement solution has 
planar components given by Eq. (A. 6) and a vertical component given by Eq. (A.7) 
while the rotation vector has planar components given by Eq. (A.9) and a vertical 
component given by Eq. (A.7) (Sokolnikoff (1946) and Knopp (1996)). 

A. 3 The three 2D examples of rectilinear defects 

In this section we consider two typical multivalued analytic functions log{Z) and 
Zlog{Z). Starting from the general uniform strain expressions Eq. (A. 5) or Eq. (A.8) 
it is easily observed that any of the holomorphic functions /" (with single- 
valued), g" and h' can provide a solution to the 2D problem. Since these functions 
can be expanded in Laurent series: 

/"(Z)=^a„Z" , ,g"(Z)=^6„Z" , /i'(Z)=^c„Z", 

— OO — OC — oo 

inside their respective convergence annuli, primitivation shows that 

+ 00 

/(^)= E f _^^, Z"+^-ci-2HZ)+a-iZ\n{Z)+A^Z + Ao, 

^ (n + l)(n + 2) 



(n + l)(n + 2)' 



h 

5(^)= E ( I IV ,^. Z"+^~b-2HZ) + b-iZln(Z) + B^Z + Bo, 



— OC 
+ 00 

E ^^"+'+C-iln(Z) + Co, 
^ — ' n + 1 

— oc 

with a_i real in order that be uniform. The relevant cases are those which give 

rise to a dislocation or a disclination, i.e. such that the functions /,/',<?', or 
are multivalued. Hence, in order to obtain non-vanishing rotation or displace- 
ment jumps, one needs to consider the following cases: 

f{Z) = -a_2ln(|) + a_iZln(f ), O-i € R, (A.IO) 
g{Z) = &-i^ln(|), (A.U) 
h{Z) = c_iln(f), c_ieiR, (A.12) 
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where i? is a constant length and to which any purely elastic term may always be 
added. In fact, from Eqs. (A. 6), (A. 7), (A. 7) and (A. 9), and from the definition 
B* := B* + iB*y with Bl given by Eq. (3.1f ), it follows that: 

B* : = Bl- iB* = [ul] - i[uy] - Ql{iz) _ 

= ^[f]-^{zif] + [g']} + Wmf'}i 

It should immediately be noted that Q,* vanishes identically since h' cannot be mul- 
tivalued. ^From Eqs. (A.10)-(A.12) and some easy computations, the only possible 
solutions are given by the following proposition. 

Proposition A.l For a straight defect line L in 2D elasticity, there are no more than 
three distinct defect classes. The two dislocation classes are the screw dislocation with 
a vertical Burgers vectors Bl, as generated by the analytical function h (f = g ~ Q), 
and the edge dislocation with a planar complex Burgers vector Bl + iBy, as generated 
by the analytical function g (f = h — 0). There is a single class of disclinations, 
the wedge disclination, which has a vertical Frank vector S7* and is generated by the 
analytical function f (g — h — 0). These functions are: 

E* ^* Z 

WEDGE DISCLINATION f(Z) = -z\n( — ) 

■'^ ' Stt ^ R' 

E*{Bl+iBl) ^, ,Z, 

EDGE DISLOGATION g(Z) = ^ -Zln( — ) 

2(1 + i>*)-K R 
iEB*-. , 

SGREW DISLOCATION h(Z) = ^ ln( — ). 

^ ' 2ir{l + v) ^R' 

For the edge dislocation, a detailed derivation is given by Eshelby (1966). 
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